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Abstract 

We investigate the entanglement between individual field theory modes in 
finite-density systems of interacting relativistic and non-relativistic fermions in 
one spatial dimension. We calculate the entanglement entropy for a single field 
theory mode and the mutual information between any two modes. The calcu- 
lation is perturbative in the four-fermion (two-body) coupling, with the leading 
contribution at order A^log(A^). At this leading order, the perturbative ex- 
pression for the entanglement entropy of a mode diverges logarithmically as the 
momentum of the mode approaches the Fermi surface from above or below. The 
mutual information between modes is largest for pairs of modes just above and 
below the Fermi momentum. The entanglement properties of modes near the 
Fermi surface are qualitatively the same if the field theory is cut off to eliminate 
modes away from the Fermi surface. 



1 Introduction 



The physics of finite density systems of interacting fermions plays a crucial role in our 
understanding of a wide variety of condensed matter systems, from ordinary materials 
to nuclear matter in a neutron star. To understand even crude macroscopic properties 
of these systems, quantum mechanics (e.g. the Pauli Exclusion Principle) is essential. 
Nevertheless, our theoretical investigations often focus on classical observables such as 
thermodynamic quantities, correlation functions, and response functions, since these 
quantities are simpler to access via experiment. In this paper, we instead introduce and 
investigate some intrinsically quantum observables in simple examples of finite-density 
fermion systems. 

One of the key features that distinguishes quantum systems from classical systems is 
the possibility of entanglement between different degrees of freedom. This entanglement 
can be quantified: given any subset A of degrees of freedom, the von Neumann entropy 
Sa = — tr(pAlogpA) of the density matrix pA = tr^ |\E')(\E'| provides a measure of 
the entanglement between A and the rest of the system A. This is known as the 
entanglement entropy of the subsystem A (for a review, see [1]). Entanglement entropy 
and related observables have been studied extensively in quantum field theory and 
many-body systems over the past several years (see, for example [21 El El S]), but these 
studies typically choose A to be the subset of degrees of freedom inside a particular 
spatial region. In this work, we study the entanglement entropy for even simpler 
subsystems: we take A a subset of momentum space, focusing on the simplest possible 
subset consisting of a single field theory mode (i.e. a single allowed momentum). 

Momentum-space entanglement entropy in quantum field theory systems was in- 
vestigated in detail recently by Balasubramanian and two of the present authors in 
[6]. There, it was emphasized that entanglement entropy in momentum space vanishes 
in the ground state of non-interacting systems and remains finite in the continuum 
limit, in contrast to the position-space entanglement entropy which is non-zero for 
non- interacting systems and diverges in the continuum limit. The work [6] showed 
further that in the presence of weak interactions, the momentum-space entanglement 
entropy can often be computed in perturbation theory, as we review in section 2. The 
goal of this paper is to carry out such perturbative calculations of the entanglement 
entropy for the simplest possible interacting finite-density fermion systems, a gas of free 
non-relativistic or relativistic fermions in one dimension perturbed by a four-fermion 
(two-body) interaction. 

We focus on two quantities in particular: the single-mode entanglement entropy 
S{p), and the mutual information I{p,q) = S{p) + S{q) — S{p,q) that measures en- 
tanglement and correlations between two individual modes at distinct momenta. In 
the non-relativistic case, for either lattice fermions (section 2) or continuum fermions 
(section 3), we show that the leading perturbative expression for S{p) diverges loga- 
rithmically in \pF — p\ as p approaches the Fermi momentum from above or below. The 
logarithmic divergence also indicates a breakdown of perturbation theory when the 
momentum is within e~^^^ of the Fermi momentum, where c is some A-independent 
constant. For continuum fermions, we also calculate the leading perturbative contri- 
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bution to the mutual information between any two modes, both in the non-relativistic 
case and (in section 5) for Dirac fermions with a {ipip)"^ interaction. This quantity shows 
discontinuities when either of the momenta cross the Fermi surface and is largest when 
both momenta are near the Fermi point. 

The renormalization group picture of such interacting fermion systems (see [7] for 
a review) suggests that the low-energy physics of the systems we consider should be 
described by a scale-invariant Luttinger liquid system. Luttinger liquids correspond to 
stable RG fixed-points, so the low-energy physics should be largely insensitive to the 
details of physics for modes far from the Fermi surface. As a check of this, we show 
that the behavior of the entanglement entropy for modes near the Fermi-surface is the 
same in a theory with a cutoff — kp] < A. As discussed in |6], there may be a 
direct connection between the behavior of systems under renormalization group flows 
and the momentum- dependence of entanglement observable; investigating this further 
is an interesting question for future work. 



2 Momentum-space entanglement in perturbation 
theory 

In this section, we review the general result [B] for entanglement entropy at leading 
order in perturbation theory and its application to the calculation of various measures 
of entanglement for subsystems corresponding to subsets of modes in momentum space. 



General result for entanglement entropy in perturbation theory 

Consider any quantum system with Hilbert space "H = "Ha ®'Hb and Hamiltonian 
H = Ha ® 1 + 1 (g) Hb + ^Hab- Letting \n) and |A^) be energy eigenstates of Ha and 
Hb respectively, a state \n) ® |A^) (an energy eigenstate of the A = Hamiltonian, 
e.g. the vacuum state) has no entanglement between the subsystems. Turning on 
the interaction, the perturbed eigenstate may be calculated using ordinary quantum- 
mechanical perturbation theory. From this, we can compute the density matrix for 
subsystem A and the associated entanglement entropy Sa- As shown in (jSj), the 
leading order perturbative expression in the non-degenerate case is 

^. = -A^log(A^) K^^^'^-^^'"'^)'^ +0(A^). (1) 

„,_^„^;vYAf ^^"^ + En - En' - En'Y 

This involves a sum over matrix elements of the interaction Hamiltonian between the 
original state and states for which both subsystems A and B have been changed. 



Entanglement entropy for a region of momentum space in quantum field 
theory 

For quantum field theory, we can start at finite volume so that the the Hilbert space 
has a discrete Fock-space decomposition 1-i = ^pHp and the unperturbed Hamiltonian 
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is a sum of terms 

p 

each acting on a single factor of the tensor product. Here, a labels the species of 
particle if there is more than one, and Ep^a is the energy of a particle of species a with 
momentum p. This may include the contribution of a chemical potential added to give 
a ground state with finite density 

m = ll<^c.,\o). (2) 

i 

We can take A to be some subset of the allowed momenta for the theory, i.e. all 
modes with a particular set of allowed wavelengths, and consider the entanglement 
entropy of the modes in region A for the state ([2]). Taking the large volume limit 
with the region A of momentum space fixed, we find [6J that the formula ([T]) gives an 
entanglement entropy for modes in the region A that is extensive (i.e. proportional to 
spatial volume), with 

Here, the matrix element Ai ji is defined by 

{^f\Hj\^,) = {2nY6{pf - Vi)Mf,{vi. ■■■.Vn) (4) 

where |\E'/) are occupation number basis elements of the form aj^^ a^' ' ' '^p„ "^^^ 
sum and integral are over the possible states appearing in the matrix element. 
Specifically, the sum is over the possible number of particles of each type present in the 
state \^ f)i while the integral is over the momenta Pa of the particles that have been 
added/removed from the initial state to produce 1^1//), with the constraint that 
we have added/removed at least one particle with momentum in the region A and at 
least one particle with momentum in the complementary region A. 

Single-mode entanglement entropy in quantum field theory 

In the case where A corresponds to a single mode, the entanglement entropy is finite 
and volume-independent in the large volume limit (since we are no longer keeping the 
momentum-space volume of the the region A fixed in the limit). We find 

S(p) = -A=log(A^)^ y"'n|^(2.mj>/-J>.)(|^ + 0(A^). (5) 

where now in the sum/integrals over the basis state |\E'/) the constraint is that the 
occupation number of the mode with momentum p is different than for and the 
occupation number of at least one mode with some other momentum has changed. 
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For the special case of spinless fermions, the Hilbert space associated with a single 
mode is only two-dimensional, and this allows us to give an explicit result for the order 

terms in the entanglement entropy in terms of the order A^log(A^) piece. From 
equation (18) in [B], we see that for 

Sip) = -\Hog{\^)a + OiX^), 

we must have 

S{p) = -A^ log(A2)a + X^a{l - log(a)) + 0{X^), 
Below, we will write explicitly only the leading (9(A^log(A^)) terms. 

Mutual information between modes in quantum field theory 

Finally, we will consider the mutual information between two modes with momentum 
p and q. Letting Ap and Ag represent subsets of momentum space corresponding to 
infinitesimal volumes d'^p and d'^q about momenta p and q, we find that the mutual 
information I{Ap, Ag) = S{Ap) + S{Ag) — S{Ap U Ag) is proportional to volume in the 
large volume limit, and also proportional to d'^p and d'^q. If we define X(p, q) by 

^(^.'^.)/^=(|^(05^fe^)' (6) 

then we find 

X(p,g) = -A^log(A^)5^ f n ^A^-r^iPf - P^) ^E^l%^ + ^(^') ' 

where now the state |\E'/) in the matrix element is required to differ from the state |\E'j) 
in the occupation numbers of modes p and q and at least one other mode. 

3 Entanglement entropy for lattice fermions with 
nearest neighbor interactions 

We now study the entanglement between modes in several models of fermions at finite 
density, starting with a systems of spinless fermions on a lattice in one spatial dimension 
with nearest neighbor interactions. 
We choose a Hamiltoniaiil] 

^Here, we have rescaled H to be dimensionless. To restore the physical dimensions, we can multiply 
by l/(a^m), where a is the lattice spacing and m is the effective particle mass, so that excitations 
with momentum p <^l/a about the unfilled state have energy p^/(2rn). 
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for which the ground state has half-fihing and an exact particle-hole symmetry [7j. We 
define the momentum-space modes ipp by 



Er a dp 
J-IL 2,71 



ijpa 



such that 

{V'p,V'J} = (27r) S2.{p-q) , 

a delta function with period 27i/a, where a is the lattice spacing. 
The Hamiltonian becomes 

dp A dp 



where 



" dPdQdpdq6{P + Q-P- g)V^l,V^J,^p^,e'(^-'')" . (9) 



The third term in ([8]) represents an adjustment to the chemical potential such that the 
state remains at half-filling in the presence of the interaction. 

For A = 0, the mode energy — cos (pa) is negative for \p\ < 7r/(2a), so the ground 
state is 

1^.) = n ^» ■ 

We can now calculate the entanglement entropy of a single mode with momentum k/a 
using ([5]). In this case, all states {"^ f) for which the matrix element (jlj) is nonzero 
have the same number of particles as {"^i), with two particles removed inside the Fermi 
surface (at momenta p and q), and two particles added outside the Fermi surface (at 
momenta P and Q). Denoting such a state by |\E'j; P, Q,p, q), we have 

(^i; P, Q, p, = -a{27r)6{P + Q-p- g)(e*(^-'')" - e*^^"^)'^ - e^^^"^)" + e'^^-P^"") 

so that the squared matrix element appearing in is 

\MiP,Q,p,q)\^ = a^lie'"^" - e'^'')(e''" - e'^")]' . 

When the momentum k/a is inside the Fermi surface, we can take p = k/a and then 
integrate over all q inside the Fermi surface and (P, Q) outside the Fermi surface. 
Performing the Q integral using the delta function sets Q = k/a + q — P in the 
integrand and the remaining integral is over all possible q inside the Fermi surface and 
P outside the Fermi surface such that Q = k/a + q — P is also outside the Fermi surface. 
The result is that the leading perturbative contribution to the entanglement entropy 
is 

c r,^ 1 x2i r^2^ f ..r, 4(cos(P - A:) - cos(P - g))^ 

^m(fc) = --^A^log(A^) / dqdP- 



477^ Jji^ (cos(g) + cos(/c) — cos(P) — cos{k + q — P)Y ' 



5 



Figure 1: Leading perturbative contribution to single mode entanglement entropy S{k) 
as a function of mode momentum for lattice fermions at half filling. This diverges 
logarithmically at the Fermi points k = ±7r/2. 



where the region of integration is 

Ri = {-k < q < --,-7, + ^ < P < --}U {-- < q < -k; - < P < ^ + tt} , 

and we have absorbed a factor of a into the integration variables. 

Similarly, when the momentum k/a is outside the Fermi surface, we can take P = 
k/a, and integrate over Q outside the Fermi surface and [p, q) inside the Fermi surface. 
The result is 



1 \2i /^2^ f 4(cOs(p-fc) -COs(p-Q)) 



Soutik) = -^X'log{X') [ dQdp 
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R^ (cos((5) + cos(A;) - cos(p) - cos(A; + Q - p)y ' 



where 



R2 = {^<Q<n;--<p<^^^-n}U{-n<Q<-k;--<p<^^} 
lJ{-k<Q<--;^<P<-}. 

It is straightforward to show that Sout{k) = 5'j„(7r/2 — k); thus, the entanglement 
entropy is exactly symmetric about the Fermi surface, a consequence of particle-hole 
symmetry. 

The function S{k) is plotted in figure [H We see that this leading perturbative 
expression diverges at the Fermi momentum; the divergence is logarithmic in \k — kpl 
and remains if we include a cutoff \k — kpl < e restricting to momenta near the Fermi 
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surface. In this case, we find 



S{k) = -^A^ log(A2) In [j^^^ + C{e) + 0{k - k^), 

where C is a momentum- independent constant of order e^. 

Since the Hilbert space for a single mode is two-dimensional, the exact mode entan- 
glement entropy is bounded by log(2). Thus, the divergence in our leading perturbative 
expression indicates a breakdown in perturbation theory when the momentum is taken 
too close to the Fermi surface. Specifically, we expect that the perturbative result is 
reliable only if it is much less than one. This requires that \k — kp] ^ e^^^'^ . 

4 Entanglement entropy for continuum non-relativistic 
fermions 

We now consider the continuum limit of the previous model, obtained by taking the 
lattice spacing to zero and adjusting the chemical potential so that states up to some 
fixed momentum (independent of a) remain occupied. Restoring the overall factor of 
1 / (a^m) in the Hamiltonian, we can rewrite the interaction f|TT]) as 

TV 

Hi = -—\ dPdQdpdq6{P+Q-p-q)i:^pij}^ijpijg{e''^''-e'''''){e'''''-e'P'') . (10) 

oZTT dTfl J ._ JL 

a 

In the limit a — > 0, this gives 

Hi = -^ 1^ dPdQdpdq5{P + Q-p-q)ijUQ^P%{Q-P){(l-p) ■ (11) 

a 

Rescaling A — t- X/{am) so that the Hamiltonian is independent of a in the limit, we 
finally obtain 

H = Ho + XHi 



where 



and 



i/o = £|v.**V(|l-.) (12) 



Hi = ^^^r^ / dPdQdpdq6{P + Q-p-q)^^p^Pl^p^,{Q-P){q-p). (13) 

6ATT m J _oo 

We now follow the same steps as in the previous section to obtain results for the single 
mode entanglement entropy. For < A; < 1, we find 

siPFk) ^ -#,A^ log(A^) I r dq r dP + r dq r - 



A-kV '\J_,^J_^ ' J_, j (P-fc)2(P-g)2 
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Figure 2: Leading perturbative contribution to single mode entanglement entropy 
S{pFk) as a ftinction of mode momentum (as a fraction of the Fermi momentum) for 
weakly interacting continuum non-relativistic fermions. This diverges logarithmically 
at the Fermi points. 



For 1 < /c < 3, we find 
S{pFk) -- 



-1 .1 
dQ dp 

k J*±2 



-2-fe 



-1 



For > 3, we have 



S{pFk) 



AMog(A^) 



2-fc 



dQ 



dp 



dQ 



-2-fe 



k + Q 
2 



dp 



-1 



{p-kf{p-Qy 



{k - Q)\2p -k-QY 

{p-kY{p-QY 



In each case, we have rescaled the integration variables by a factor of to make 
them dimensionless. All of these integrals may be evaluated analytically to obtain 



where 



S{pFk) - 



Pf 



-{k + If ln(l - A;) - (1 - kf ln(l + k) + k^ + 2 ln(2)) + | + 2 ln(2) < A; < 1, 

m ^{ -(k+ ir Hk -i) + (k+ ir n + h2)) -f(k+i) + f-^^ i<k<3, 

3<k. 
(14) 

The results for /c < are obtained using S{—k) — S{k). 
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3(fc2_-l) 
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Figure 3: Regions of {k,l) space with different behaviors for X{k,l). Results for the 
unshaded regions may be obtained from the results for the shaded regions using the 
indicated symmetries of X. 

The single- mode entanglement entropy S{pFk) is plotted in figure [2] . It diverges 
logarithmically at k = ±1 with the behavior on either side described by 

f{k) = -4 In |1 - A;| + — sgn(l - k) + 41n(2) + C(l - k). 

3 

At A; = 3, the function f{k) and its first, second, and third derivatives are all continuous, 
with a discontinuity appearing only in the fourth derivative. 

4.1 Mutual information between modes 

We can also look at the entanglement structure in more detail by calculating the mutual 
information between individual modes with momenta ppk and ppl. Specifically, we 
calculate the function X{ppk,ppl) defined in ([6]). 
The mutual information X{ppk,pFl) satisfies 

X{pFk,pFl) = X{pFl,PFk) = X{-pFk, -pfI) , 

so we can restrict to the region {k > 0, |/| < k} and find X for the other values using 
the symmetries. For each choice of k and /, the integral in ([7]) is over distinct pairs of 
momenta such that together with the momenta ppk and pfI, we have two momenta 
inside the Fermi surface and two momenta outside the Fermi surface (otherwise the 
matrix element in ([7]) vanishes). Performing one integral using the delta function, we 
are left with a single integral in each case. For the various regions depicted in figure |31 
we find the following results: 

• Region A: {0 < A; < 1, |/| < k} 
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m'^ J _ 



Pf dP{k-lf{2P -k-l) 



2 



27r {P-kf{P-lY 



27rm2 1^ ' ' + U {k + l){l + l) 

Region B: {/ < -1, -/ < A; < 2 - /} 

^.PF dp{k-lf{2p-k-lf 



X{pfKpfI) = -AMog(A2)^ / 



m? Jk+i 2tt {p — k^ip — I) 
-A^MA^)TT^|2(/-A;)ln 



Region C: {-1 < / < 1, 1 < < 2 + /} 

dp{p-lf{2k-p-lf 



I{pFk,PFl) = -X'log{X')^J^ 



_i 2tt {k-pY{k-lY 
.2, .^2^ f , l{k + lf-{l + lf 



Region D: {-1 < / < 1, A; > 2 + /} 

/"^ dp{p-lf{2k-p-lf 



I{pFk,PFl) = -XHog{X')^ [ 

J _ 



2 



m2 27r {k-pY{k-l) 



27rm2 [3 (A; - /)2 P - 1 

Region E: {A; > 1, Z > 1} U {A; > 3, 2 - A; < / < -1} 

X{pFk,pFl) = . 

The function Z{pFk,pFl) is plotted in figure HI with plots for specific values of k given 
in figure |5l We see that this leading-order contribution to X is generally discontinuous 
as one momentum crosses ±Pf and diverges when both momenta approach one of the 
Fermi points, unless the two momenta are equal. 

5 Entanglement entropy for relativistic fermions 

In this section, we consider Dirac fermions with a four-fermion interaction, described 
by the action 

S = I d'^x ^^itp'j'^dfj^ip — m^jj'ip — Xip'ip'ijj'ip^ . 
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Figure 4: Mutual information X{ppk,pFl) between individual modes for non-relativistic 
fermions with Fermi momentum pp. 
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Figure 5: Mutual information I{ppk,ppl) vs / for non-relativistic fermions with k = 
(top left), k = 0.75 (top right), k = 0.95 (bottom left), and k = 1.25 (bottom right). 
The overall scale for I is arbitrary 
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This reduces exactly to the non-relativistic model in the previous section (with A = A/2) 
if we take chemical potential n = m+fi with jl <^m and consider observables related to 
energy scales small compared to m so that the particle number is fixed and antiparticles 
decouple. 

Here, we consider general values of the chemical potential /i and the corresponding 
ground state for which all particle states with energy less than /i are occupied (and 
none of the antiparticle states are occupied). 

The field may be expanded as usual in terms of creation and annihilation operators 

as 

z^(x) = j (a, uu e-^^- + hi v, e^'^) . (15) 

where and hk respectively annihilate a particle and antiparticle with momentum k. 

For the calculation of any of the entanglement observables in section 2, the non- 
vanishing matrix elements are those between the ground state and states obtained 
by 

1. adding two particles and two antiparticles, corresponding to the operator combi- 
nations a^h'^a^h'^ in the expansion of J dx {ipipY; 

2. adding two particles and removing two particles, corresponding to a^aa^a] or 

3. adding two particles, removing one particle and adding an antiparticle, corre- 
sponding to a^aa^l)^ and a^Wa^a. 

Particles can only be added outside the Fermi surface and can only be removed inside 
the surface; antiparticles can be added anywhere, but cannot be removed without 
annihilating the ground state. 

Using ( JTSjl . we can write the relevant terms in the interaction Hamiltonian in terms 
of the momentum-space creation and annihilation operators as 

^2 = 2 y I J] ^^'^^ 1 w(A;i)t;(A;2)w(A;3)^(fc4)4i 4.2 43^^4 (27r)(^(A:i + + k^, - k^) 

^ y jn 2n^%j' '} ^(^1)^(^2)^(^3)^(^4)41 43Qfc2Qfc4(27r)(5(A;i + fcs - ^2 - k^j 

Using these, we can calculate the matrix elements appearing in the calculation of 
entanglement observables. As an example, consider the matrix element of Hi between 
the ground state and a state \fJ',PiP2',P3P4:) where two particles with momenta pi 
and P2 and two antiparticles with momenta p^ and p4 have been added to the ground 
state. We find 

{lI,PlP^;P2Pi\Hi\lj) = ^^^^^^2{u{pi)v{p2)u{p3)v{p4) - U{pi)v{p4)u{p3)v{p2)) 

(16) 
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so that from (jl]), 



we get 



1 \u{pi)v{p2)u{Pz)v{Pi) - U{pi)v{pi)u{p'i)v{p2)\' 

{pi ■ p3 - m^){p2 ■ P4 - m^) 



UJiUJ2UJzUJi{uJi + CU2 + 1^3 + 

Similarly, we can calculate the matrix element of H2 between the ground state 
and the state \fi,PiP3P2',P4:) where two particles and an antiparticle have been added 
with momenta pi, ps and p^ respectively, and a particle with momentum p2 has been 
removed. We find 



\M 



P1PW2\PA I 



1 



\u{pi)v{pi)u{p'i)u{p2) - U{ps)v{p4)u{pi)u{p2)\' 



4UIUJ2UJ3UJ4 



+ m 



^1^21^31^4 (1^1 + ^3 - ^2 + 
J2iPhP2,P3,P4) ■ 



Finally, we have 



\u{pi)u{p4)u{ps)u{p2) - U{ps)u{p4)u{pi)u{p2)\' 



{AEf 



4:UiU2UJ3U4 



{pi ■ ps - m^){p2 ■ Pi - nr) 

UOiUJ2^Z^a{0Ji + - 072 - 1^4) 
^3(Pl,P2,P3,P4) 



{Ui + CJ3 - W2 - ^^4)' 



for the matrix element of between the ground state and the state |/i,PiP3P2P4) where 
particles have been added with momenta pi and pa and particles with momenta p2 and 
Pi have been removed. 



Single-mode entanglement 

To calculate the entanglement entropy for a single mode with momentum p, we use 
the expression ([5]), taking the sum over the three types of final states discussed above. 
However, we find that the integrals in the terms involving Hi and H2 diverge. Thus, as 
noted in [B] for the case without chemical potential, the leading perturbative expression 
for the single-mode entanglement entropy in this model is ill-defined. Since the exact 
answer is necessarily less than log(2), this must indicate a breakdown of perturbation 
theory, as discussed in more detail in section 4.4 of [6]. Thus, for this model, we focus 
on the mutual information between modes, which can be computed in perturbation 
theory. 



Mutual information between modes 

To calculate the mutual information between modes, we use ([7]). In each case, X{p,q) 
is calculated using matrix elements for which the occupation number of the particle 
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modes with momenta p and q have been changed relative to the ground state and for 
which occupation numbers for two other particles or antiparticles (with momenta P 
and Q) have been changed. 

When p and q are both inside the Fermi surface, we have: 



X(p, q) 



-A^log(A^)- 



\jjP j dQ6{P + Q-p-q)MP,p,Q,q) 



where and indicates that the integration variable ranges outside and inside the 
Fermi surface, respectively. 

When p is inside, and q is outside: 



X(p, q) 



Ztt 



j dP j dQ 5iP + Q - p + q)J2{q,p, P,Q) 
+ JdpJ dQ6{P-Q-p + q)J3{q,p,P,Q) 



This also covers the case when p is outside and q is inside, since I{p, q) = I{q,p). 
Finally, when both p and q are outside: 



X(p, q) 



dP j dQ5iP + Q + p + q)Jiip,P,q,Q) 

j dP j dQSi-P + Q+p + q)J2ip,P,q,Q) 

dP [ dQ5i-P-Q + p + q)J3{p,P,q,Q) 
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All these integrals are straightforward to perform numerically. We find behaviour 
qualitatively similar to the non-relativistic case, with discontinuities at the Fermi mo- 
menta and the largest mutual information when both modes are close to the Fermi 
point. As an example, in figure E] the mutual information X{p,q) at /i = 10m is plot- 
ted as a function of q for several fixed values of p. These may be compared with the 
non-relativistic results in figure [51 

For comparison, we plot in figure [7] the mutual information X(p, q) vs q for several 
values of p in the case with zero chemical potential, where the unperturbed ground 
state is the Fock-space vacuum. Here, we find smooth behavior with the largest mutual 
information between pairs of momenta of opposite sign and momenta of order the mass 
scale. The mutual information falls off as 1/g for fixed p. 
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Figure 6: Mutual information X(p, q) vs q for relativistic fermions with /i — 10m and 
p — (top left), p — 7.5 (top right), p — 9.5 (bottom left), and p — 12.5 (bottom 
right). Momenta are given in units of m, so the Fermi points are at ±10. The overall 
scale for I is arbitrary. 
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Figure 7: Mutual information X(p, q) as a function of p and q for relativistic fermions 
with = 0. Momenta are given in units of m. 
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